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Abstract. A stationary magnetic dipole immersed in an electric field carries 
"hidden" mechanical momentum. However, the fate of this momentum if the fields 
are turned off is unclear. We consider a charge-and-dipole hidden momentum 
configuration, and turn off the fields by collapsing a null shell onto the system, 
forming a black hole. In numerical calculations we find that the black hole receives 
a kick corresponding to 0.1% of the initial stored momentum. When extrapolated to 
apply to purely gravitational phenomena, this efficiency suggests a role for the hidden 
momentum kick mechanism in generating the binary black hole "superkicks" observed 
in numerical simulations of Einstein's equation. 
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1. Background and Motivation 

First described over fifty years ago pQ , the phenomenon of hidden mechanical momentum 
continues to generate interest [21 EJ SI El EJ [T] . The basic idea is simple. In flat spacetime, 
any time-independent, conserved stress-energy T^ u must have zero total momentum P % = 
J T l0 d 3 x (when this integral is defined). For matter coupled to electromagnetic fields, we 
then have for stationary configurations that the total mechanical plus electromagnetic 
field momentum is zero, 

PM + pEM = Q (1) 

(Here the total stress energy T^ u = T*5 + T^ U M is given by a matter component Tj^ 
and the field component T^J* , and the momenta P M and P EM are defined in terms of 
their respective stress tensors.) However, using stationary sources it is easy to create 
stationary field configurations whose field momentum is non-zero. Thus, the sources 
contain non-zero momentum — there is mechanical momentum in bodies at rest! This is 
the hidden momentum. 

The original and simplest example is a magnetic dipole immersed in an electric 
field. Consider a point dipole of strength jl located at a position Xdi P in a stationary 
electric field E = — V$ assumed to vanish at infinity. Then the field momentum is 

P em = ±-JexB = -±- I V$xfi=i-| rtxS = J*f=-jlx£\^. (2) 

The first equality is the definition of P EM ( the second uses E = — V$, the third 
integrates by parts, the fourth uses the Maxwell equation V x B = AirJ, and the 
fifth uses the point-dipole current J = —jl x VS(x — Xdi P ). By equation 0, the hidden 
mechanical momentum is then 



5Ai 



tfx^kp- (3) 



Equation ^ gives the mechanical momentum present in a stationary dipole immersed 
in an electromagnetic field. How is this momentum stored? For a simple model of a 
current loop dipole, the electric field modifies the flow of current in such a way that 
the (relativistic) particle momenta of the charge carriers sums precisely to the simple 
formula above [21 IE] • 

The above discussion applies to stationary situations only, and it is natural to 
ask the fate of the hidden momentum if the fields are turned off or the configuration 
is otherwise destroyed by a time-dependent process. Is some (or all) of the hidden 
mechanical momentum converted into ordinary momentum, resulting in a kick (balanced 
by emitted radiation)? Or does the system just relax to the trivial configuration 
containing a body at rest? The process is difficult to model, since one must include 
the agent turning off the fields as part of the calculation, and the result will in general 
depend on the details of the turning off process. To our knowledge, no such calculation 
has yet been carried out. A recent study of hidden momentum [3] concluded that while 
some impulse may be delivered, "there is no obvious reason why this impulse should 
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equal the momentum originally stored in the fields — all we can say in general is that 
the total momentum afterward, like the total momentum before, is zero." 

In this paper we will determine the fraction of momentum that is converted into 
impulse in one particular scenario — collapse to a black hole. Specifically, we will consider 
a point electric charge placed near a point magnetic dipole, both surrounded by a shell 
of matter assumed large and light enough that its gravity does not (initially) distort 
the electromagnetic field configuration. Thus, initially, we have a perfectly ordinary 
hidden momentum configuration. However, at a time to the shell begins to collapse 
at the speed of light, beginning the process of black hole formation. By the "no hair 
theorem" the final black hole can only support a simple monopolar electric field, so the 
final configuration has no hidden momentum and no field momentum. By monitoring 
the momentum radiated (relative to the rest frame of the black hole) we determine 
that a kick corresponding to 0.1% of the initial stored momentum is delivered to the 
black hole. To perform these calculations we adapt a model due to [9], which treats 
the electromagnetic fields as small perturbations to the collapse spacetime and exploits 
causality to simplify the calculations. Although the "form a black hole" method of 
turning off the fields may seem unduly complicated (or completely out of the blue), this 
type of trick makes it likely to be one of the simplest situations to model. 

A completely separate motivation for studying black hole formation in the presence 
of hidden momentum comes from numerical simulations of black hole mergers. A 
surprising discovery [101 [111 EE] was that for certain configurations of black hole binaries 
containing spinning black holes, there is a bobbing up and down of the entire orbital 
plane that correlates to a large kick velocity obtained at merger. In previous work 
[I], this behavior was explained using the concept of hidden momentum. By revealing 
the bobbing to be associated with a kinematical effect of spin that is present even for 
systems for which there can be no kick (an example is spinning balls connected by a 
string), the authors of [4| concluded that the kick should not be viewed as an inertial 
continuation of the bobbing. But what, then, causes the kick? It was noted that, at least 
in an electromagnetic analog of orbiting charged magnetic dipoles, there is an orbitally 
modulated field momentum present outside the bodies, along with corresponding hidden 
momentum within the bodies. The authors suggested that a rapid merger event allows 
the release of some of this field momentum, leading to a kick that will depend sinusoidally 
on the phase at merger, as observed in numerical simulations. While this picture 
accounts qualitatively for the observed behavior, there was no quantitative estimate 
of the percentage of the field momentum that a merger could release. Furthermore, 
there was no means of distinguishing the proposal from a competing scenario — also 
qualitatively successful — in which the kick is caused by the bobbing motion's modulation 
of the dominant orbital emission [13] . The calculations of this paper provide an estimate 
of the efficiency of the momentum release mechanism in a context where the competing 
mechanism is absent, thereby helping to distinguish the two scenarios. The implications 
for our understanding of black hole merger kicks of are discussed in section [3j 
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2. Model and Results 

In their early investigation of the fate of perturbing 
fields during gravitational collapse, De La Cruz, 
Chase, and Israel [H] introduced the null shell 
collapse spacetime pictured in figure [2j The 
spacetime is constructed by joining flat spacetime 
(region I) to Schwarzschild spacetime (regions 
II, III, IV) on a spherical hypersurface (blue line) 
that is initially tangent to the direction of time 
symmetry (i.e., stationary), but at a later time 
suddenly becomes tangent to the ingoing radial 
null direction (collapsing at the speed of light). 
Following the standard junction conditions [TU 
ITS] , the induced metric is required to match on 
the shell, and the jump in extrinsic curvature is 
interpreted as the stress-energy of the shell. This 
spacetime models a sudden and rapid gravitational 
collapse of a shell of matter. In the words of the 
authors of [9], the model is "highly artificial from an astrophysicist's point of view, but 
does not violate any of the principles of relativity theory." 

We work in coordinates (t, r, 8, <p) with d t the direction of time symmetry, r the 
areal coordinate for symmetry two-spheres, and (6, <p) spherical coordinates on those 
spheres. For t < t , the metric is that of the well-studied static spherical shell, given in 
our coordinates by 




Figure 1. The collapse spacetime. 
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dr 2 + r 2 dtl 2 , 



r < r 
r > r , 



(4) 



where dfl 2 = d6 2 + sin 2 9d(fi 2 is the metric of the two-sphere. Here ro > 2M is the initial 
coordinate radius of the shell. However, we will take ro to be very large, 



r > M, 



(5) 



so that for t < to the spacetime is simply Minkowski spacetime everywhere. 

On this Minkowski spacetime we place a point charge q and a point magnetic dipole 
fi. The magnetic dipole is located at the origin and is oriented in the z (8 = 0) direction, 
while the charge is located on the x axis (<ft = 0) at a distance D < 2M from the origin. 
(We require D < 2M so that the charge ends up inside the black hole.) In a suitable 
gauge the associated fields are 



e=o m= 



V A-7T — 



l> 



sin 2 9, 



(6) 
(7) 
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where we have expanded the point charge field A t in spherical harmonics Y £m (with 
the conventions of [20J). Here r < (r > ) denotes the lesser (greater) of r and D, and an 
overbar denotes complex conjugation. By assumption, equations ^ and give the 
fields for t < to- However, by causality they must also give the fields in the entirety 
of regions I and II, and, by continuity of the electromagnetic field, and can therefore 
provide characteristic]!] initial data for an evolution into the remaining portions of the 
spacetime. Since our main concern is with the momentum radiated to infinity, we confine 
the evolution to region III, the Schwarzschild exterior. We use double-null coordinates, 
related to Schwarzschild coordinates (t, r) by 

u = t - r* (8) 

v = t + n, (9) 

where r* = r + 2Mlog(r/(2M) — 1) is the tortoise coordinate. The two-sphere at the 
onset of collapse (meeting point of regions I, II, III) is described by u = uq = to — r*(ro) 
and v = v Q = t Q + r*(r ). For t > t , v = v is the collapsing shell's surface, while 
u = uq represents a spherical shock front propagating outward at the speed of light, 
bringing news of the sudden onset of collapse. The (future) event horizon is located at 
u — > oo (fix v), while v — > oo (fix u) represents future null infinity the vantage-point 
of a distant observer receiving the radiation. The shell, shock front, horizon, and null 
infinity bound our evolution region III, u < u < oo and v < v < oo. The initial data 
is given on the shell and shock front, v = vq and u = u . 



2.1. Field Equations 

To solve the time-dependent Maxwell equation on the Schwarzschild background, we will 
use the master function approach [T6J, [T7j as formulated in [18] . The vector potential is 
expanded into four scalar amplitudes, A^ 71 , v lm and v £m , 

oo I 

A a = K m Y em (10) 

1=0 m=-i 
oo I 

a a = J2Y1 { vimz T + v im x e ^) , (ii) 

i=\ m =~e 

where indices a,b,c, ... run over the coordinates t,r and indices A, B,C... run over the 
coordinates 9, </>. Here Z^ 1 = OaY 1 " 1 is the even parity vector spherical harmonic, while 
X 1 ™ = eA B dBY £m is the odd parity vector spherical harmonic. (The Levi-Civita tensor 
€ab is that of the unit two-sphere (e^ = sin 9), and the index is raised with the metric of 
the unit two-sphere.) One then defines master functions ip± in terms of the amplitudes 

| Characteristic initial refers to data on the union of two intersecting null surfaces, in this case the 
surfaces u = uq and v = vq. 
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by 

^ = r 2 (d t A e r m - d r A £ t m ) (12) 
^ = £(£ + l)v £m , (13) 

and, using Maxwell's equations, it turns out that these functions entirely determine the 
field strength by 

00 1 \.lm. 

£=0 m=-£ 
oo £ 1 

fu = EE 777TTT i^- mx A m + (1 - ^)^ m 4 m } (15) 

£=1 m=-£ l{l + ^ 

oo £ 

f r i = EE l^- m ^ m + (1 - (16) 

1=1 m=-£ 1 ^ + ^ 

oo e 

F ^ = J2J2 ^- m sin^ m . (17) 

1=1 m =-£ 

Back in terms of double-null coordinates, Maxwell's equations imply that the master 
functions satisfy 

d 2 %b± ( 2M\£(£ + 1) . . 



Equation (18) applies only to the modes with £ > 0; the £ = mode is non-radiative and 
requires a separate treatment. (The non-radiative character follows immediately from 
Gauss' law. Since our interest is in the radiation, we do not consider this modej§]) Note 



that equation (18) is the same for all m and each choice of + or — . The initial data 
for our evolution is provided by the static solutions, equations @-(j7])- Since D < 2M, 
the initial data surfaces u = uq and v — v o satisfy r > D, and we require the point 
charge solution only for r > D (i.e., r> = r, r< = D). Computing the master functions 
associated with the initial data then gives 

^ = 4nq^— - Y im {\ , 0) (19) 



2£ + 1 \r 

. •w^S if£ = landm = 
= <J V 3 r ' _ (20) 

0, otherwise 
Here we have used the fact that = a/3/(47t) sin 2 6. 

§ Our initial solution for £ = 0, A t oc 1/r, is also a solution in Schwarzschild spacetime (when areal 
coordinates are identified, as they arc in our shell metric). This observation allows one to see explicitly 
that the I = mode does not evolve when the shell begins to collapse. 
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Equations (19)-(20) provide initial data on u = uq and v — Vq for the evolution of 
equation (18). In order to avoid redundant integrations over different m- values (which 
satisfy the same wave equation), it is convenient to introduce real amplitudes f± such 
that 

Att — 

^ = 4^^(1,0) (21) 

V>i° = /I^§*. (22) 
The amplitudes satisfy the same wave equation, 

d 2 fi ( 2M\£(£ + 1), 



+ -4lT^/± = ° (^>«). ( 23 ) 



dudv \ r J 4r 
with initial values given by 

fl\u=u ,v=v = q(£ + l) (jj (24) 

/_ \u=u ,v=v = ~- (25) 

2.2. Radiated Energy and Momentum 

Our main interest in this evolution is the total energy and (especially) momentum 
radiated away. At future null infinity, v — > oo, the spacetime is flat and we may 
compute the total flux of energy and momentum using expressions valid in Minkowski 
spacetime. Let t a and (xi) a (i = 1, 2, 3) be the time and space translation Killing fields 
of flat spacetime. (In Cartesian Minkowski coordinates, we have t a = (1, 0, 0, 0) and 
(xiY = S^i.) These vector fields, when dotted into a conserved stress-energy tensor, 
give the associated current of energy and momentum, which can be integrated over the 
surface v — > oo to give the total radiated energy and (i th component of) momentum. 
For the total electromagnetic energy and momentum lost through future null infinity in 
our spacetime, we have 

AE = lim / (\T a Bt p d a v) r 2 dttdu (26) 

v — >oo J 

AP { = - lim I (\T a B{xifd a v) r 2 dQdu, (27) 

v— >oo J 

where 6, <fi range over the two-sphere (area element dfl = sin OdOdifi) and u runs from — oo 
to oo. (However, for our field configuration the flux vanishes for u < Uq and one may 
begin the integral at u = u .) Here is the stress-energy tensor of the electromagnetic 
field, 

47T7V = F m F v a - \g^F aP F^. (28) 

The radiated energy takes an elegant form in terms of the master functions; using 
the orthogonality of the spherical harmonics it is straightforward to show that 

AE = E £ ^Tiy ™ J (l ^ m|2 + '^- m|2) (29) 
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where the bars denote complex norm. In deriving this equation and analogous flux 
equations below, we have used the fact that as v — > oo the general solution of equation 



(23) is the sum of a pure function of u and a pure function of v. In terms of our real 



amplitudes, we have 




The general expression for the radiated momentum (analogous to equation (29) 
for the energy) will be less elegant and more laborious to derive, since it will involve 
Clebsch-Gordon-type coupling coefficients. (Each spatial Killing vector (xi) a is i = 1, 
so that triple products of spherical harmonics appear under the integral.) However, 
from the symmetry of the problem and the "addition of angular momentum" rules, it is 
straightforward to establish the following: (a) the radiated momentum in the x direction 
takes the form of a coupling between even parity modes with £- values separated by 1; (b) 
the radiated momentum in the y direction will come entirely from the I = 1 even-parity 
mode multiplied by the £ = 1 odd-parity mode; and (c) the radiated momentum in the 
z direction is vanishing. Evaluating the angular integrals explicitly for low values of £, 
we have 

AP x = -\im J (^(dJl)(d u fl) + ^(d u f 2 + )(d u fl) + ...)du (31) 

AP y = - lim / (l(d u f l _)(d u fl)) du (32) 
AP Z = 0. (33) 

The radiated momentum in the x direction would persist in the absence of the magnetic 
dipole and can be thought of as resulting from the asphericity of the collapse of a shell 
onto an off-center point charge. This asphericity will result in a slight asphericity of 
the emitted radiation, leading to a kick. On the other hand, the radiation in the y 
direction requires both the dipole and the charge and may be thought of as due to 
a release of the initial field momentum (entirely in the y direction) during the rapid 
collapse event. Since the initial data for each real amplitude f± scales as (D/r) e , we 
can anticipate that this will be the dominant contribution, i.e., \P y \ > \P X \- (We find 
that \P y \ m lQqD 2 / (fiM)\P x \.) Furthermore, it is clear that we lose little accuracy by 



dropping terms higher order in £, represented by the ... in equation (31). (We find 
that the second term is smaller than the first by a factor of about 75 for order unity 
parameters.) It is clear that we may therefore neglect the third and higher terms, 



already not displayed in equation (31). 



2.3. Parameter Scalings 

Some further simplification is helpful for numerical purposes and for revealing the scaling 
of the results with the parameters of the model. Since f} and f]_ are proportional 
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initially (see equations (24) and (25)) and satisfy the same linear wave equation (23), 
they are in fact proportional for all time. We may thus keep track of both 1 = 1 
amplitudes by a single fiducial amplitude /. We will define / to be the numerical 
function obtained by integrating the wave equation (23) with i = 1 and M = 1 and 
with initial data given by / = 1/r (on u = Uq and v = Vq in the limit of a large shell, 
u — > — oo). Using equations (23)-(25), the momentum flux formula (32) becomes 



AR 



where we define the fiducial integral 7 by 



lim 



' M 3 ' 



(d u f) 2 du, 



(34) 



(35) 



with the factor 2/3 introduced for convenience. The scaling of AP y with /i, q, and D is 
evident from equations (24) and (25), while the factor of M -3 could have been guessed 
on dimensional grounds|j|] The integral 7 is simply a number that may be computed 



numerically. From equation (35) we have that 7 is manifestly positive, so that the kick 
is guaranteed to be in the direction of the initial hidden momentum, as expected. We 
find numerically that 7 ~ 1.35 x 10~ 3 (see section 2.4). By conservation of momentum, 



1.35 x 10 3 (see section 
the net kick delivered to the final black hole is 



.kick 



APy/M 



fiqD 



7. 



(36) 



From equations (|2|)-(j3|) , the initial field configuration has momentum P y M = fJ-q/D 2 
(counterbalanced by hidden mechanical momentum p™ ech = —fiq/D 2 ). We may define 
the efficiency e^ck of the hidden momentum kick as the fraction of this field momentum 



that is radiated away during the collapse. From equation (34) we have 



Ckick 



pEM 

y 



D 

M 



I. 



(37) 



Thus for separation distances D of order the mass M of the shell, the efficiency is about 
0.1%. For the Maximum allowed separation distance, D = 2M, the efficiency is 1.0%. 

If we keep only the t = 1 terms in the energy flux, then equation (30) shows that 
the fiducial integral 7 also controls the energy flux, 



A75I 



[i 2 + q 2 D< 
M 3 



7. 



(38) 



We find numerically that the I = 2 term is a correction of fractional order 10~ 3 (D/M), 
while the I = 3 term is a correction of fractional order 10 _5 (7J/A7). Thus equation (38) 
should give an excellent approximation to the Energy flux. 



|| To establish equation (34) it is helpful to introduce coordinates u = u/M and v = v/M, which agree 
with u and v in the fiducial case M = 1. 
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Finally note that we may divide equations (34) and (38) to obtain 



AEl- 



fiqD 



fi 2 + q 2 D 2 ' 



(39) 



The fiducial integral has canceled out, and we have a fully analytic result for the ratio 
of the ?/-momentum flux to the energy flux. We see that the ratio is controlled by 
a symmetric combination of the initial magnetic dipole moment \i and electric dipole 



moment qD, reaching a maximum of 1/2 when the moments are equal. Equation (39) 



could also have been derived directly from equations (30) and (31) using the fact that 
/| = {2qD/n)fl. 

To summarize, we first presented our problem in terms of a general master-function 
formulation of Maxwell's equations on a Schwarzschild background. Then we reduced 
that prescription to one involving real amplitudes /£, whose initial values are given in 



equations (24) and (25). Once equation (23) is integrated with these initial values, the 



fluxes of energy and momentum are given by equations (30)- (33). Finally, we introduced 
a fiducial /, corresponding to integrating equation (23) with M = 1 and initial data 1/r. 



This integration should be carried out with initial shell radius taken to be very large, 
u — > — oo. The integral /, defined by equation (35) in terms of this /, then contains all 



the information about the dominant energy flux AE\i=\ (equation (38)) and momentum 



flux AP y (equation (34)) for any choice of model parameters q,fi,D,M. 



2.4- Numerical Implementation 



We evolve equation (23) using a 2nd order finite difference discretization, 



f N = fE + fw-fs- h 2 V c - (f E + f w ) 



(40) 



1 on a grid of 
1, 2, 3. In principle 



where the indices N, S, W, E refer to north, south, west, east on the grid as in reference 
[2"T] . h is the grid spacing, and V c is approximated by V c ~ (1/2) (V^ + Vw) and 
V = (1 - 2M/r)£(i + l)/(4r 2 ). We evolve equation Q with M 
u G [Mo,M max ] and v G [0, v max \ with initial data l/r e for £ 
our model requires uq = — oo, w max = +oo and t> max = +oo, and we settle on values 
Uq = —50, M max = 150 and v max = 200 as large enough that all observables have 
converged. Figure [2] shows the result of these evolutions with / evaluated on f max , 
approximating future null infinity J? . We performed each simulation for three different 
resolutions with grid spacing of h = 0.1, 0.2, 0.4. The dashed line shows the convergence 
factor c = log 2 ((f^h — /2/1) / (/2ft — fh)), which should be c = 2 for 2nd order convergence. 
Notice that the waves reach null infinity around u = —10, peak around u = +10, and 
decay by about u = 60. Since our shell is at v = 0, these retarded times correspond 
(respectively) to emission from r* = {5,-5,-30}, or r w 4.53M (initial emission), 
r w 2.05M (peak of emission) and r w 2.00M (end of emission). Thus the vast majority 
of the emission comes from the strong-field regime. 
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Figure 2. The field f(u, v max ) as a function of u on the grid boundary v = v max 
which approximates J? . The field (left scale) is evolved for the three different cases 
I = 1,2,3. The dashed line shows the convergence factor c s» 2 (right scale) and 
demonstrates good 2nd order convergence. 



The y momentum flux and the dominant energy flux are both captured by the 



fiducial integral I, equation (35), which is computed from the i = 1 evolution. We find 

J=1.35xl(r 3 . (41) 

In establishing this result we have quantified three sources of error: finite resolution 
(e res = 10~ 7 ), finite extraction radius (e Vmax = 10 -5 ), and finite initial shell radius 
(e Uo = 3 x 1CT 5 ). Our resolution error estimate is done via Richardson extrapolation. For 
the finite extraction radius error, we vary t> max and fit a function of the form Io + h/v max 
in 1/iVax- F° r the initial shell radius error we similarly fit a function of the form I +Ii/uq 
in 1/uq. The value I quoted above is extrapolated in v max , but not in resolution h or Uq. 
From the i = 2 and i = 3 evolutions we also compute the higher order corrections to 



the radiated energy, equation (30), and the radiated x momentum, equation (31). For 
order unity parameters, the subleading energy term is a fractional correction of order 
10 -2 , while the x momentum flux is smaller than the y momentum flux by a factor of 



order 10. The text below equation (33) gives the precise numbers, along with parameter 
scalings. 

As a further check of our code we reproduce the results of the numerical run done 



Hidden Momentum and Black Hole Kicks 



12 



in reference j9]. In particular, for their parameter choices of M — 1/2 and uq = —10, we 
are able to reproduce their integral Ii (see text below their equation (2)) to the accuracy 
given in the paper. 

3. Discussion 

The main result of the paper is the y kick velocity of the black hole, 

|< c 1^ (400 km/s) (42) 

where we now measure the kick in kilometers per second (while [i,q,D,M are still in 
Gaussian geometrized units). This result holds precisely within our model, which we 
emphasize is a well-defined physical setup that can in principle be realized in nature. 
Since it is unlikely to be realized in practice, however, its main utility is in revealing 
the character and size of the basic physical effects it entails — in our case, the idea that 
a large kick can be produced when a quasi-static configuration with hidden momentum 
is suddenly destroyed. While one would expect on basic physical grounds that some 
kick is produced, we see no a priori argument for the size, and indeed our investigation 
could in principle have revealed the effect to be negligible. Our numerical calculation of 
the prefactor of 400 km/s establishes that the hidden momentum mechanism can be a 
factor in the generation of large kicks. 

3.1. Application to Gravity 

In order to compare directly to the kicks observed in binary black hole simulations we 
must make several extrapolations. First, we must use the analogy between gravitation 
and electromagnetism, involving q ~ m and // ~ S, where m and S are the mass and 
spin of a black hole in the (comparable mass, comparable spin) binary. Second, since the 
shell of mass M enforces our merger and produces a final black hole of mass ~ M, we 
must also take M ~ m to approximate the merger of black holes of mass ~ m. Third, 
to approximate the premerger configuration of black holes, we should take the distance 
between the particles to also be of order the mass/size of the black holes, D ~ m. Finally, 



we must extrapolate the formula (42) beyond the test-Maxwell-field regime where it is 
strictly valid, since the above parameter identifications entail q ~ M and fi ~ M 2 . The 
combination of these extrapolations — in which we have left out all numerical factors of 
order unity — may migrate the kick velocity by an order of magnitude. If so, then we 
predict a kick in the range of 40 to 4000 km/s, the upper end of which does correspond 
to the observed superkicks. 

We therefore believe that the hidden momentum kick mechanism may have a role to 
play in explaining the superkick phenomenon, and could possibly account for the entire 
effect. However, an alternative proposal p3] also predicts kicks of the correct qualitative 
and roughly the correct quantitative character. In this scenario one imagines that the 
bobbing motion of the black holes provides a small modulation to the dominant emission 
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(due to the quasi-circular motion), which induces an asymmetry in that emission, 
resulting in the kick. Since the bobbing motion is of direct relevance in producing 
the kick in this scenario, we will call this the "motion scenario". This scenario does 
not rely on the presence of a merger to generate the kick, and a number of simulations 
of grazing collisions of black holes (i.e., near the threshold between capture/merging 
and scattering) [22] have indicated that large kicks can be obtained even in the non- 
merging case. However, in these simulations the black holes complete at least a full orbit 
before separating [23], and one could therefore still view the kick as a consequence of 
release of quasi-static field momentum due to the sudden separation of the bodies. Thus 
the simulations significantly inform the discussion but do not conclusively distinguish 
between the two scenarios. 

The work reported here helps further distinguish the scenarios by isolating the 
relevant physics of our scenario (a rapid event destroying the quasi-static configuration) 
and creating a situation where the relevant physics of the motion scenario is entirely 
absent. Indeed, our "binary" consists of particles at rest, with the emission entirely 
determined by the merger's destruction of the initially static field configuration. That 
we can recover kicks in the ballpark of those seen in numerical simulations without 
any orbital motion at all suggests that the rapid destruction of the quasi-static field 
configuration may indeed be playing an important role in producing the superkicks. 
However, the uncertainty inherent in the extrapolation of our results to the gravitational 
case leaves plenty of room for the effect to be too small to fully account for the superkick. 
The motion effect could be dominant, or both effects could be playing an important role. 

3.2. Hidden Momentum More Generally 

We conclude with a discussion of other physical and astrophysical systems in which 
hidden momentum may have a role to play. First, since the phenomenon is basically 
independent of the internal structure of the body (being enforced by the presence of 
quasi-static field momentum and the conservation of total momentum), similar kick 
phenomena are to be expected for mergers and near-merger scattering of neutron stars, 
boson stars, or any other object that could be simulated. However, the kick will be 
smaller for bodies that are less compact. Returning to the electromagnetic case where 
the notion of hidden momentum is most precisely defined, one could imagine that 
a magnetized neutron star moves through an external magnetic field, generating an 
"external" electric field in the local rest frame. This would be an astrophysical example 
of the original hidden momentum setup, and our heuristic would predict a kick if the 
neutron star underwent collapse. However, even assuming the strongest known magnetar 
moving through the field of another such object at an orbital separation of ten times 
the radius, the total hidden momentum corresponds to a kick of only about a millimeter 
per second. (Furthermore, plasma dynamics would likely be important for this system.) 
Having a large electromagnetic hidden momentum kick requires compact objects with 
strong intrinsic electric and magnetic fields. One possible avenue for further exploration 
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would be the addition of electric charge to the spinning black holes already considered, 
which generates an effective magnetic dipole moment as well. In this case there should 
be significant hidden momentum of both the electromagnetic and gravitational variety. 
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